Lattice-Independent Approach to Thermal Phase Mixing 
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We show how to achieve lattice-spacing independent results in numerical simulations of finite- 
temperature stochastic scalar field theories. We generalize the previous approach of Ref. Jl| by 
obtaining results which are independent of the renormalization scale. As an application of our 
method, we examine thermal phase mixing in the context of Ginzburg-Landau models with short- 
range interactions. In particular, we obtain the lattice-spacing and renormalization-scale indepen- 
dent critical value of the control parameter which determines the free-energy barrier between the 
two low-temperature phases. We also propose a simple procedure to extract the critical value of 
control parameters for different choices of lattice spacing. 



PACS numbers: 05.10.-a; 05.70. Jk; ll.10.Wx 
I. INTRODUCTION 

During the past decade, the study of equilibrium and 
nonequilibrium dynamics of field theories has greatly 
benefited from the widespread availability of worksta- 
tions capable of millions of floating point operations per 
CPU-second. 

One of the most popular applications of computers in 
the physical sciences is the examination of phenomena 
which are generated by nonperturbative effects. These 
include nonlinear dynamical systems with a few, several, 
or an infinite number of degrees of freedom. Of these, 
we are particularly interested here in the latter, as they 
represent a unique challenge to computational physics. 
Implementing field theories in the computer implies dis- 
cretizing not only time but also space: the system is cast 
on a finite lattice with a discrete spatial step, effectively 
cutting off the theory both in the infrared (by the lat- 
tice size) and in the ultraviolet (by the lattice spacing). 
Although in classical field theories an ultraviolet cutoff 
solves the Rayleigh- Jeans ultraviolet catastrophe, the so- 
lution comes with a high price tag: whenever there is 
dynamical mixing of short and long wavelength modes, 
the results will in general depend on the shortest distance 
scale in the simulation, the lattice spacing. To be sure, 
in many instances this dependence on small spatial scales 
does not affect qualitatively the physics one is interested 
in: for example, very near criticality for Ising systems, 
where spatial correlations in the order parameter diverge 
[]2J, or are controllable in some way |3|. However, in many 
other cases one is interested in achieving a proper con- 
tinuum limit on the lattice which is independent of the 
choice of ultraviolet cutoff. These include a wide range 
of phenomena which have triggered much recent interest, 
from pattern formation in fluid dynamics ||] to simula- 
tions of phase transitions and topological defect forma- 
tion, which often use stochastic methods pj. 

In the present work, we are concerned with curing, 
or at least greatly alleviating, the lattice-spacing depen- 
dence that appears in stochastic simulations of scalar 
field theories. We will show that it is indeed possible 



to obtain results which are lattice-spacing independent, 
as long as proper counterterms are added to the lat- 
tice effective potential. Following a suggestion by Parisi 
|6| , lattice-spacing independent results were recently ob- 
tained by J. Borrill and one of us within the context 
of finite-temperature symmetry restoration in a simple 
Ginzburg-Landau model 0. However, that study fo- 
cused on a regime where the large temperatures needed 
for symmetry restoration compromised the approach to 
obtain lattice-spacing independent results, which is based 
on a perturbative expansion in powers of the temper- 
ature. Furthermore, no attempt was made to obtain 
results which were independent of the renormalization 
scale. Thus, in that study, the numerical prediction for 
the critical temperature depends on the particular choice 
of renormalization scale. 

Here, we would like to apply an expanded version of the 
method proposed in Ref. Q to a related problem, phase 
mixing in Ginzburg-Landau models. The distinction be- 
tween phase mixing and symmetry restoration is made 
clearer through the following argument. Suppose a sys- 
tem described by a Ginzburg-Landau free energy density 
with odd powers of the order parameter <p(t, x) is rapidly 
cooled from a high temperature to a temperature where 
two phases can coexist. The system was cooled so as to 
remain entirely in one of the two phases. The odd term(s) 
could be due to an external magnetic field (linear term), 
or to the integration of other fields coupled to <j>, as in 
certain gauge theories (cubic term) M , or in de Gennes- 
Landau models of the nematic-isotropic transition in liq- 
uid crystals ||. Due to the odd terms, there is a free- 
energy barrier for large-amplitude fluctuations between 
the two phases. [Of course, small-amplitude fluctuations 
within each phase are also possible, but less interesting.] 
This barrier is usually controlled by the coefficients of 
the odd terms in the Ginzburg-Landau model, which we 
will call the control parameter(s). 

Suppose now that the system is held at the temper- 
ature where the two phases have the same free energy 
densities (sometimes called the critical temperature in 
the context of discontinuous phase transitions) and that 



we are free to change the value of the control param- 
eter (s). The question we would like to address is how 
does the system behave as a function of the control pa- 
rameter^), that is, as the free-energy barrier for large- 
amplitude fluctuations between the two phases is varied. 
As is well known, the mean-field theory approach breaks 
down when large-amplitude fluctuations about equilib- 
rium become large enough. Thus, we should expect that 
the prediction from the Ginzburg-Landau model, that 
the system remains localized in one phase until the bar- 
rier disappears (when the control parameter goes to zero) 
will eventually be wrong. There will be a critical value for 
the control parameter beyond which nonperturbative ef- 
fects lead to the mixing of the two phases. [Note that due 
to the odd terms, there is no symmetry to be restored.] 
In the language of the Ginzburg-Landau model, the sys- 
tem should at this point be described as having a single 
well, centered at the mean value of the order parameter. 
We would like to obtain the lattice-independent critical 
value of the control parameter for this phase mixing to 
occur. 

It is important to distinguish between phase coexis- 
tence and phase mixing. As is well-known, phase co- 
existence will generally occur when a system is cooled 
into the so-called phase coexistence region of the phase 
diagram. In this case, the system will relax into its low- 
est free-energy configuration via spinodal decomposition. 
Here, we are preparing the system initially outside the 
phase coexistence region, namely in one particular phase 
only. In the infinite- volume limit, mean field theory pre- 
dicts the system will remain there, since, as the two min- 
ima are degenerate, the nucleation of a critical droplet 
would cost an infinite amount of free energy. Phase mix- 
ing is a nonperturbative phenomenon characterized by 
large-amplitude fluctuations not included in the mean- 
field approach. It signals the breakdown of mean-field 
theory. 

This paper is organized as follows: In the next section 
we describe the continuum model we use and some of its 
properties. In section 3 we describe the lattice implemen- 
tation and how simulations using a bare lattice potential 
give results which depend severely on the lattice spac- 
ing. In section 4 we show how to cure this dependence, 
and also how to make the simulations independent of the 
choice of renormalization scale. We conclude in section 
5 with a brief summary of our results and a discussion of 
future work. 



II. THE MODEL IN THE CONTINUUM 

Our starting point is the 2-dimensional Hamiltonian, 
(we use c = fcs = 1) 



where the homogeneous part of the free energy density is 
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This choice of V{4>) is inspired by several models of nu- 
cleation in the condensed matter H and high-energy 
physics literature, in particular in recent models of the 
electroweak phase transition M. The several parame- 
ters in V((j>) allow one to apply it to several situations 
of interest. However, we note that here the order pa- 
rameter is a scalar quantity, and thus the critical be- 
havior of this model belongs to the universality class of 
the 2-dimensional Ising model ]10|,p|. It is quite straight 
forward to generalize our results to systems in different 
numbers of spatial dimensions. 

At the critical temperature T c 2 = T 2 2 /(l - 2a 2 /9a\), 
the system exhibits two degenerate free energy minima 
at 



= and 
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while at the temperature T 2 the barrier between the two 
phases disappears. Throughout this work, we will be in- 
terested in the behavior of the system at T c . One reason 
for this choice has to do with the use of a perturbative 
expansion which is in powers of T; at T c the expansion pa- 
rameter is sufficiently small, allowing us to stay at 1-loop. 
Another reason is that we are interested in measuring the 
breakdown of mean-field theory in terms of parameters 
controlling the free-energy barrier, and the calculations 
are much simpler at T c , as we will see next. 

According to the model described by Eq. 0, at T c , 
unless a = (or A — > oo) there will always be a barrier 
separating the two phases: at T c , this model does not 
predict phase mixing to occur. It is thus very convenient 
to introduce the shifted field 
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and write the shifted homogeneous free energy density as 
(dropping the primes) 



V Q {4>) = - 7 V(T)</> 2 + -0 4 + A(T)</> + constants , (5) 
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The shifted free energy density is just the usual Ginzburg- 
Landau free energy density with an external magnetic 
field A(T). Note that A(T C ) = and the two minima 



are degenerate, as they should be. We now introduce 

the dimensionless variables 9 = T/T2, t = \faTit, x = 

y/aT 2 x, 4> = -$jr, according to which we can write, at 
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Since we will keep the system at 9 c (a, A), the only two 
control parameters are a and A. In what follows, we will 
fix A = 0.1 for simplicity. This was also the choice in a 
previous study of phase mixing in the same system, which 
did not address the issue of lattice-spacing dependence 
p[ . We will also drop all tildes, except in the plots 
and captions, where unshifted, dimensionless variables 
are used and marked explicitly. 



III. NUMERICAL RESULTS: BARE LATTICE 

A. Description of the Simulation 

As mentioned in the introduction, we would like to 
study the behavior of the system described in the previ- 
ous section when coupling to an external thermal bath 
promotes fluctuations about equilibrium. We will con- 
sider the situation where the system is initially prepared 
in the phase given by <fi — in the unshifted potential or, 
more generically, the left well. Since we are only inter- 
ested in the final equilibrium value of the system, we will 
simulate the coupling of the scalar field <f> to the thermal 
bath using a generalized Langevin equation, 



d 2 4> 2 50 dV« 



(10) 



where the viscosity coefficient 77, set equal to unity in 
all simulations, is related to the stochastic force of zero 
mean £(x, t) by the fluctuation-dissipation relation, 



(£(x, i)£(x',t')) = 2t]95{x - x')6(t - t') . 



(11) 



The system is discretized and put on a square lattice 
with side length, L, equal to 64 for all the simulations, 
but several lattice spacings, Sx, and time steps, 5t, are 
used. For Sx =1.0, 0.8, and 0.2 the respective time steps 
are St = 0.2, 0.1, and 0.02. We have, of course, checked 
the stability of the program for these choices of lattice 
parameters. Using a standard second-order staggered 
leapfrog method (which is second order in both space 
and time) we can write, 



vm+ 2 l + ^St 

<f>i,m+l = (f>i,m + St<j) im+ i 



(12) 



where i-indices are spatial and m-indices temporal, over- 
dots represent derivatives with respect to t and primes 
with respect to <f>. The discretized fluctuation-dissipation 
relation now reads 



it- /- \ 0/1 ®i,j Om n 



so that 
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(13) 



(14) 



5x 2 St l ' m 

where Gi^ m is taken from a zero-mean unit-variance 
Gaussian. 



B. Results from Bare Lattice Simulations 

Keeping the system always at the critical tempera- 
ture 6 C , we are interested in its behavior as the free- 
energy barrier between the two equilibrium phases is 
changed. We will measure the value of the ensemble- 
averaged and area-averaged order parameter ((j))A{t) = 
1/A J d 2 xcf)(x, t) for several choices of the lattice spacing 
Sx, taking note of its final equilibrium value, <p cq . In fig- 
ure 1 we show the results for ((/>)^(t) for several choices 
of lattice spacing and a = 0.45. The dependence on 
lattice-spacing is quite evident; different lattices produce 
different physics. 
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FIG. 1. ((/>)a(£) for the bare potential, for several choices 
of lattice spacing, for a = 0.45. 
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FIG. 2. Phase diagram for the bare potential for various 
lattice spacings. 

In figure 2, we show the phase diagram depicting phase 
mixing as a function of a for different choices of the lat- 
tice spacing Sx. The phase diagram is constructed by 
defining the "phase-mixing order parameter" , 



5<f,(a) = |0 O q - 0max|/</>ii 



(15) 



where 



LX = a9 c /3\ is the location of the max- 
imum of the free energy density separating the two 
phases. Clearly, as a decreases, the free-energy barrier 
decreases and larger-amplitude fluctuations between the 
two phases become more probable. Below a critical value 
q c , cj> ccl just tracks the location of the maximum, indicat- 
ing complete phase mixing, or the breakdown of the mean 
field theory of Eq. |[ 

The problem, though, is that phase mixing, or the 
breakdown of mean-field theory, occurs for values of a c 
which are strongly dependent on the value of Sx, as can 
be seen from figure 2. For the range of 5x investigated, 
0.2 < Sx < 1, we obtained 0.355 < a c < 0.40. In the 
next section, we argue that this dependence can be ef- 
fectively cured by including proper counterterms to the 
lattice potential. 



IV. APPROACHING THE CONTINUUM ON 
THE LATTICE 

A. Computing the Lattice Effective Potential 

Setting up a continuum system on a lattice introduces 
two artificial length scales, the ultraviolet momentum 



cutoff A = tt/Sx and the infrared momentum cutoff 
Al = tt/L, where L is the lattice size. In the contin- 
uum limit, L — > oo, and Sx — > or, equivalently, the 
number of degrees of freedom TV = (L/Sx) d — > oo. The 
coupling to the thermal bath induces fluctuations at all 
allowed length scales. We should thus expect that the 
lattice simulation is related to a continuum model with 
both infrared and ultraviolet cutoffs. In order to obtain 
the lattice effective potential, we start by analyzing the 
divergences of the related continuous model. 

For classical field theories in 2 dimensions, the corre- 
sponding 1-loop corrected effective potential is given by 



Vil(0)=W) + ~ 
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where the primes denote derivatives with respect to <j>. 
Performing the integration and making all variables di- 
mensionless we obtain, 



Kl(0) 



(17) 

The infrared cutoff does not introduce a divergence as 
Al — > 0, but it does introduce finite corrections to 
Vil, or finite size effects, which become small as L in- 
creases. These become more severe near criticality, but 
well-known scaling behavior can be used to regulate this 
dependence |12j| . As we will further argue below, for our 
purposes we can safely set A^ = 0. This is not the case 
for the ultraviolet cutoff. The reader can see now why it is 
useful to use the shifted potential of Eq. || as opposed to 
the original one of Eq. Et all divergences are quadratic in 
4>, simplifying the computations considerably, while the 
physical results, of course, remain unchanged. This is 
why we added only the counterterm Be/) 2 above. 

The counterterm B is computed by imposing the renor- 
malization condition 

V(' L (cf> BN )=V£(cf> BN )=M 2 , (18) 

where M is the arbitrary renormalization scale and we 
[Note that M here is dimension- 
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less (tilde is dropped), being defined as M — M/T2] One 
obtains, 
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(19) 
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Applying this to the shifted potential of Eq. |L we 
obtain, for the 1-loop renormalized continuum potential, 
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Recall that at 9 C the linear term proportional to A{&) 
vanishes. Since the counterterm cancels the dependence 
on the ultraviolet cutoff, we define the lattice effective 
potential as jj| 



Vi at t(^) = Vo{cf>)+B(M)cj ) 2 



(21) 



In figure 3 we show the results of repeating the simu- 
lations of figure 1 but now adding the counterterm to 
the lattice simulations following Eq. |l[ The addition of 
the counterterm practically eliminates the lattice-spacing 
dependence of the results. Figure 3 also shows the near 
elimination of lattice-spacing dependence for a = 0.40. 
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FIG. 3. (</>}a(?) with the counterterm added and M — 1, 
for several choices of lattice spacing, for a = 0.45 and a — 0.40 



B. Extracting the Critical Value of the Order 
Parameter 
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FIG. 4. Phase diagram for M = 1. 
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FIG. 5. Phase diagram for M = 10. 



In figures 4 and 5, we show the phase diagrams using 5</, 
defined in Eq. na as a function of a for different choices 
of lattice spacing 8x. These are to be compared with 
figure 2. Figure 4 is for a choice of renormalization scale 
M = 1, while figure 5 is for M = 10. It is clear that the 
results for different lattice spacings converge around one 
value of a c . 



We compute a c as follows: for a given value of a we 
perform several (« m ax) measurements of cf) eq by varying 
the lattice spacing, which we call <^'(a). Their aver- 
age is simply (4> cq (a)) = y^ max 0* q (a) /i max , while 
the departure from the average for each measurement 

is, A</4j = \4>i q - (4> eq }\/(4> cq ). 



Near criticality, the results are naturally poorer due 
to the existence of long-range correlations in the field. 
We can use this fact to our advantage, since we expect 
that, at criticality, the departure from the average defined 
above is maximized, that is, the quantity 



<A0 eq (a)> 



Er a^. 



cq 



(22) 



reaches a maximum at a c . This can be clearly seen from 
figure 6 for the same choices of lattice spacings (or coarse- 
graining scales) as in figures 4 and 5. 
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FIG. 6. (A</> eq (a)} for M = 1 and M = 10, with the 
respective & c 's at the maxima. 

The measured value of a c is now a c ~ 0.365 ± 0.005, 
for M = 1, and a c ~ 0.435 ± 0.005 for M = 10. 

We have thus achieved lattice-spacing independence on 
the measurement of a c . Clearly, the error in a c could be 
further decreased by taking a larger number of measure- 
ments of 



-cq- 



However, since our main goal here is to 
show the convergence of the results for different lattice 
spacings, we are not concerned with very high-accuracy 
measurements. Nevertheless, the values for a c still de- 
pend on the renormalization scale, which is arbitrary. In 
the next subsection, we show how to obtain lattice results 
which arc independent of M . 



C. Achieving independence of renormalization scale 
on the lattice 

As with conventional renormalization theory, the 
renormalized potential should not depend on the choice of 



renormalization scale |13|] , One usually solves the renor- 
malization group equations to find how the couplings 
vary with the scale. Here, we propose a simpler approach 
which works quite well on the lattice implementation of 
scalar field theories. It is an interesting question how to 
generalize it to more complex models. 

Consider the 1-loop renormalized potential Vjjf (<j)) as 
given in Eq. Efl. The superscript M is a reminder that 
this potential is renormalized at a given scale M. Now 
consider an equivalent potential renormalized at another 
scale M', Vfjf (<f>). Since the divergences are quadratic, 
this potential has a shifted mass /i' 2 . By imposing that 
the two potentials are identical, V^((j>) = V^ ((/)), we 
obtain a condition on the shifted mass // , approximating 
In (-// 2 + 3AM') ~ In (-/i 2 + 3AM) , 
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Thus, we can always relate a theory with a choice of M 
to any other theory with M' by redefining the mass /s 2 
according to Eq. gJJ We claim that this is also the case 
for the lattice effective potential. 

As an illustration, we show the phase diagram for 
M' = 10 in figure 7, where the results for M' = 10 
were obtained after scaling fj 2 according to Eq. B3 in 
the lattice potential of Eq. ElL It is practically indistin- 
guishable from the phase diagram for M = 1 shown in 
figure 4. 
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FIG. 7. Phase diagram for M' = 10. 



Figure 8 demonstrates clearly that M' = 10 has the 
identical a c previously found for M = 1, within our level 
of accuracy. This is in stark contrast to figure 6, where 
the values of a c for M — 1 and M = 10 were very differ- 
ent, as evidenced by its "twin peaks" structure. 
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rameter around criticality for different choices of lattice 
spacing (or coarse-graining scales). Thus, we effectively 
turn a weakness of lattice simulations into a strength, 
something that can be useful for the examination of crit- 
ical phenomena of continuous field theories in fairly small 
lattices. 

We plan to expand the present study to investigate 
the effects of spatio-temporal memory on the dynamics of 
nonequilibrium fields. Recent results have shown that the 
effective Langevin equation for self-coupled scalar sys- 
tems exhibits colored and multiplicative noise Q. It is 
possible to expand the two-point function characterizing 
the noise (or noises) in terms of a "persistence factor", 
which defines short or long-term memory, spatial, tem- 
poral, or both. The possible impact of this kind of noise 
on the nonequilibrium dynamics of fields remains largely 
unexplored. 
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FIG. 8. (A0 eq (fi)) for M = 1 and M' = 10 showing the 
same value of a c . 



V. SUMMARY AND OUTLOOK 

We have investigated the continuum limit of lattice 
simulations of stochastic scalar field theories. In par- 
ticular, we have proposed a method to obtain not only 
lattice-spacing independent results, but also results inde- 
pendent of the renormalization scale of the lattice effec- 
tive potential. We illustrated our approach by examining 
a Ginzburg-Landau model which exhibits phase mixing 
depending on the values of the parameters controlling 
the free-energy barrier for large-amplitude fluctuations 
between the two low-temperature phases in our model. 
Thermal fluctuations of the order parameter are induced 
by coupling it to a thermal bath at fixed temperature T c , 
defined as the temperature where the two phases have 
the same free energy density. We simulate the dynam- 
ics using a generalized Langevin equation with Gaussian 
noise, which brings the system to its final equilibrium 
state. 

The results were presented in terms of phase diagrams 
which clearly illustrate the effectiveness of our approach. 
We also proposed a simple way of determining the critical 
value of the control parameter for phase mixing, which 
uses the spread in values of the equilibrium order pa- 
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